Phonon renormalisation in doped bilayer graphene 
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We report phonon renormalisation in bilayer graphene as a function of doping. The Raman G 
peak stiffens and sharpens for both electron and hole doping, as a result of the non-adiabatic Kohn 
anomaly at the Y point. The bilayer has two conduction and valence subbands, with splitting 
dependent on the interlayer coupling. This results in a change of slope in the variation of G peak 
position with doping, which allows a direct measurement of the interlayer coupling strength. 
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Graphene is the latest carbon allotrope to be discov- 
ered [IJ [2] |31 H] [5] . Near-ballistic transport at room tem- 
perature and high carrier mobilities [2J [3J |U [5J [5J |7J [H], 
make it a potential material for nanoelectronics [3 ITUl 
11 , especially for high frequency applications. It is 
now possible to produce areas exceeding thousands of 
square microns by means of micro-mechanical cleavage of 
graphite. An ongoing effort is being devoted to large scale 
deposition and growth on different substrates of choice. 

Unlike single layer graphene (SLG), where electrons 
disperse linearly as massless Dirac fermions[U[21[3l[4l[5], 
bilayer graphene (BLG) has two conduction and valence 
bands, separated by 71, the interlayer coupling [T!Zl I13j. 
This was measured to be~0.39eV by angle resolved pho- 
toelectron spectroscopy [14] ■ A gap between valence and 
conduction bands could be opened and tuned by an ex- 
ternal electric field (^lOOmeV for^l0 13 cm~ 2 doping)[15 
H5] , making BLG a tunable-gap semiconductor. 

Graphene can be identified in terms of number and 
orientation of layers by means of clastic and inelas- 
tic light scattering, such as Raman [17] and Rayleigh 
spectroscopies [T51 H5]. Raman spectroscopy also allows 




FIG. 1: (color online). Experimental setup. The black dot- 
ted box on Si02 indicates the polymer electrolyte (PEO + 
LiC104). The left inset shows an SEM image of the SLG and 
BLG. Scale bar: 4pm. The right inset, the 2D Raman band 



monitoring of doping and defects [H 120"! I2T1 1221 1251 I2H 125] . 

Indeed, Raman spectroscopy is a fast and non-destructive 
characterization method for carbons [26 . They show 
common features in the 800-2000 cm -1 region: the 
G and D peaks, around 1580 and 1350 cm -1 , respec- 
tively. The G peak corresponds to the E 2g phonon at 
the Brillouin zone center (r). The D peak is due to 
the breathing modes of sp 2 atoms and requires a defect 
for its activation [271 123 HH]- The most prominent fea- 
ture in SLG is the second order of the D peak: the 2D 
pcak[17j. This lies at ~ 2700 cm -1 and involves phonons 
at K+Aq[T7]|23]. Aq depends on the excitation energy, 
due to double-resonance, and the linear dispersion of the 
phonons around K[T71 [251 HO]- 2D is a single peak in 
SLG, whereas it splits in four in BLG, reflecting the evo- 
lution of the band structure [T7J. The 2D peak is always 
seen, even when no D peak is present, since no defects 
are required for overtone activation. 

In SLG, the effects of back and top gating on G- 
peak position (Pos(G)) and Full Width at Half Max- 
imum (FWHM(G)) were reported in Refs[2jJl [2TJ [25]. 
Pos(G) increases and FWHM(G) decreases for both elec- 
tron and hole doping. The G peak stiffening is due to the 
non-adiabatic removal of the Kohn-anomaly at r|20, 31 . 
FWHM(G) sharpening is due to blockage of phonon de- 
cay into electron-hole pairs due to the Pauli exclusion 
principle, when the electron-hole gap is higher than the 
phonon energy [20] [32], and saturates for a Fermi shift 
bigger than half phonon energy [SUJ, [21] |32J. A simi- 
lar behavior is observed for the LO-G - peak in metal- 
lic nanotubes|33j. for the same reasons. The conceptu- 
ally different BLG band structure is expected to renor- 
malize the phonon response to doping differently from 
SLG[T3J [33] • Here we prove this, by investigating the 
effect of doping on the BLG G and 2D peaks. The G 
peak of doped BLG was recently investigated [35]. and 
reproduced that of SLG, due to the very low doping 
range(^ 5 x 10 12 cm~ 2 ), not enough to cross the sec- 
ond BLG subband. Here we reach much higher values 
(~ 5xl0 13 cm~ 2 ), probing the further renormalisation re- 
sulting from crossing to the second BLG subband. 

We recently demonstrated a SLG top-gated by poly- 
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FIG. 2: (color online). Raman spectra of (a) SLG; (b) BLG 
at several Vtg- Red lines fits to the experimental data. 
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mer electrolyte [23] , able to span a large doping range, 
up to^5xl0 13 cm~ 2 [24 . This is possible because the 
nanometer thick Debye layer [531 E23 EZ] gives a much 
higher gate capacitance compared to the usual 300nm 
SiC>2 back gate[5]. We apply here this approach to BLG. 
Figjl] shows the scheme of our experiment. A sample is 
produced by micromechanical cleavage of graphite. This 
consists of a SLG extending to a BLG, as proven by the 
characteristic SLG and BLG 2D peaks in the inset of 
Figjljl7 . An Au electrode is then deposited by pho- 
tolithography covering both SLG and BLG, Fig[l] Top 
gating is achieved by using a solid polymer electrolyte 
consisting of LiClO.4 and polyethelyne oxide (PEO) in 
the ratio 0.12:1 [23]. The gate voltage is applied by plac- 
ing a platinum electrode in the polymer layer. Note that 
the particular shape of our sample, consisting of a BLG, 
with a protruding SLG, ensures the top gate to be effec- 
tively applied to both layers at the same time. This would 
not necessarily be the case for a monolithic BLG, where, 
due to screening effects, the gate would give a separate 
evolution of the Raman spectra of the top and bottom 
layers [55]. Measurements are done with a WITEC con- 
focal (X50 objective) spectrometer with 600 lines/mm 
grating, 514.5 nm excitation, at<lmW to avoid heat- 
ing. For a given top gate voltage, Vtg, spectra are 
recorded after 10 mins. Figs|2]a,b) plot the spectra as 
a function of Vtg- We use Voigt functions to fit the G 
peak in both SLG and BLG. The SLG 2D band is fit- 
ted to one Lorentzian. The BLG 2D band is fitted to 
four Lorentzians,2Di A ,2DiB,2D 2 A,2D 2B [17], Figjl] As 
previously discussed, two of these, 2Di^ and 2D 2 a, are 
much stronger 1 17]. Thus, we focus on these. 

To get a quantitative understanding, it is necessary 
to convert Vtg into a E^ shift. For electrolytic gat- 
ing, the chemical potential is eVrG = E F LG + ecf> SLG — 
Ef LG + ecf) BLG . The electrostatic potential <j) = ^- is 
determined by the geometrical capacitance Ctg an d car- 
rier concentration n (e is the electron charge), while Ep/e 
by the chemical (quantum) capacitance of graphene. For 



FIG. 3: (color online). Pos(G) for (a) SLG; (c) BLG as a 
function of Fermi energy. FWHM(G) of (b) SLG;(d) BLG as 
a function of Fermi energy. Solid lines:theoretical predictions. 
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SLG, n SLG = \iE F , where \x 

are spin and valley degeneracies, 7 = ^70 a, with 70 the 
nearest-neighbor tight binding parameter, a the graphene 
lattice parameter, and vp is the Fermi velocity. Thus: 



eV T G = E F + vEi 



(1) 



For BLG [HI EH SO] 



^BLG_ 



and 



,BLG 



/47i-Ef + E F ] for E F < 71 



2fj,E F for E F > 71. Thus: 



eVTG= {l + V7i)Ep 
= E F 



vE 2 F , E F < 71 
E F > 71 



(2) 



2uE*, 



where v 



We take Ctg 



TrC T G(fiv F ) 2 - " c uc ""= WJ U 2.2 X 

10~ 6 Fcm~ 2 [33 , and 7i=0. 39eV constant with doping 
(since its variation for n up to~10 13 cm" 2 is <5%~)[T31 
rTB"]). Eqs |l|2| then give Ep as a function of Vtg- 

Fig(3]plots the resulting Pos(G), FWHM(G) as a func- 
tion of E F . In SLG, Pos(G) does not increase up to 
Ep ~0.1eV (~ hujo/2), where wo is the frequency of the 
Ei g phonon in the undoped case (fiwo/(27r?ic)=Pos(Go), 
with c the speed of light), and then increases with Ep. 
Fig|3]},d indicate that in SLG and BLG, FWHM(G) de- 
creases for both electron and hole doping, as expected 
since phonons decay into real electron-hole pairs when 
E F < hu /2[20. Fig|3}; plots Pos(G) of BLG. (i) Pos(G) 
does not increase until Ep ~0.1eV (~ Hu>o/2).(u) Be- 
tween 0.1 and 0.4eV, the BLG slope R= ° F ' is smaller 
than the SLG one.(iii) A kink is observed in Fig. 3b at 
Ep ~0.4eV.(iv)Beyond Ep >0.4eV the slope is larger 
than in SLG.(v) The kink position does not significantly 
depend on 71 used to convert Vtg in Ep(e.g. a ~ 66% 
change in 71 modifies Ep by ~ 6%). 

These trends can be explained by considering the ef- 
fects doping on the phonons: (i) a change of the equilib- 



rium lattice parameter with a consequent "static" stiff- 
ening/softening, APos(G) s ';(ii) the onset of "dynamic" 
effects beyond the adiabatic Born-Oppenheimer approx- 
imation, that modify the phonon dispersion close to the 
Kohn anomalies, APos(G) dyn ^M [3T]. Thus, the total 
phonon renormalization can be written as|20| 1ST]: 

Pos(G Bf )-Pos(G ) = APos(G) = APos{G) st +APos(G) 

(3) 
For SLG, we compute APos(G) st by converting E F into 
the corresponding electron density n , then using Eq.3 
of Ref.[31 . For BLG, we assume n BLG equally dis- 
tributed on the two layers, each behaving as a SLG with 
an electron concentration n BLG /2. Eq.3 of Ref. [31] is 
then used to compute APos(G) st for BLG. APos{G) dyn 
is calculated from the phonon self-energy II|41j: 

hAPos(G) dyn = Re[IT(£; F ) - Ii{E F = 0)]. (4) 

The electron-phonon coupling (EPC) contribution to 
FWHM(G) is given by gU H21 143] : 



dyn 



FWHM(G) BPC = 2Im[U(E F ) 



(5) 



The self-energy for the E 2g mode at T in SLG is [2"Ull3"T] : 



U(E F y 



2e + Hioq + id 



|e|de, 



(6) 



while for BLG it is given by |13j : 



U(E F ) 



BLG 






s,s' j,j' 

[f( e sjk) — f(e s 'j'k)][esjk — Es'j'fc] 

(e„jfc - e s 'j'fc) 2 - (huj + id) 2 



(7) 



where a' = ^"/J^ff , A uc = 5.24 A 2 is the graphene 
unit-cell area, M is the carbon atom mass, /(e) = 
l/[exp( e r y ) + 1] is the Fermi-Dirac distribution, 5 is a 
broadening factor accounting for charge inhomogcneity, 
EPC(r) is the electron phonon coupling [46 . s = ±1 and 
,s'=±l label the conduction (+1) and valence (-1) bands, 
while j = 1,2 and j'=l,2 label the two parabolic sub- 
bands. e S jk is computed from Eq.2.8 of Ref. [13 , and <f>\-, 
is given by Eq.3. 1 of Ref. [T3]. By using Eqs 6J7in Eqs 4J5 
we get APos(G) dyn ,FWRM(G) EPC for SLG and BLG. 
To compare Eqs |3|5| with the experimental data, 
we use a' = 4.4 x 10~ 3 (obtained from the DFT 
values of EPC(r) and v F %2H [30]), the experi- 
mental huj for SLG and BLG, and T=300K. 
5 is fitted from the experimental FWHM(G) to 
FWHM(G)=FWHM(G) BPC +FWHM(G)°, with 

FWHM(G) a constant accounting for non-EPC effects 
(e.g. resolution and anharmonicity) . For SLG (BLG) 
we get 6 = 0.13eV (0.03eV) and FWHM(G) =4.3cm~ 1 
(S.lcm^ 1 ). These 5 values are then used to compute 
Pos(G). Note that the relation between n and E F 




FIG. 4: (color online). Phonon renormalization for BLG: (I) 
E F < huj , (II) huj < Ep < 71, (III) E F > 71. Blue and 
red arrows correspond respectively to positive and negative 
contributions to II. Solid and dashed arrows correspond to 
interband and intraband processes respectively. 



implies that charge inhomogeneity causes different Ej? 
broadening in SLG and BLG (e.g. <5n~10 12 cm~ 2 would 
give 0.13eV and 0.03eV in SLG and BLG, respectively). 
The solid lines in Fig. [3] are the theoretical Pos(G) 
and FWHM(G) trends. The experimental and theoret- 
ical FWHM(G) are in excellent agreement, as expected 
since the latter was fitted to the former. The theoreti- 
cal Pos(G) captures the main experimental features. In 
particular, the flat dependence for \Ef\ < 0.1 eV in both 
SLG and BLG, and the kink at - 0.4 cV in BLG. This 
kink is the most striking difference between SLG and 
BLG. It is the signature of the second subband filling in 
BLG. Indeed, a shift of E F , by acting on /(e) in Eq|7] 
modifies the type and number of transitions contributing 
to LI. The only transitions giving a positive contribu- 
tion to LT are those for which \e s j,k ~ Cs' j',k\ < ^Oi 
i.e. a subset of those between (s = — l;j = 1) and 
(s = 1; j = 1) (interband transitions, solid blue lines in 
Fig.ffl. Interband transitions with \e 8 j,k — Cs',j',k\ > ^o 
(solid red lines in Fig. El and all intraband (between 
(s = ±l;j = 1) and (s = ±l,j = 2), dashed red lines in 
Fig. El contribute to II as negative terms. It is convenient 
to distinguish three different cases: (I), \E F \ < Hojq, (II) 
huj < \E F \ < 71, and (III) \E F \ > 71. For simplicity let 
us assume E F > (the same applies for E F < 0).ln case 
(I), positive contributions from interband transitions are 
suppressed, and new negative intraband transitions are 
created. This results in strong phonon softening at low 
temperatures |35j . At T=300K, these effects are blurred 
by the fractionary occupation of the electronic states, 
resulting in an almost doping independent phonon en- 
ergy (see FigJ3J)). In case (II), a shift of E F suppresses 
negative interband contributions and creates new neg- 
ative intraband transitions. By counting their number 
and relative weight (given by <&j,ji/{e s _j t k — £ s ',j',k))> one 
can show that interband transitions outweight intraband 
ones, resulting in phonon hardening. Case (III) is simi- 
lar to (II), with the difference that the second subband 
filling suppresses negative intraband transitions at k^K, 
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FIG. 5: (a) Ratio of 2D and G peaks intensities for SLG (solid 
circles) and BLG (open circles) as a function of n. (b) Position 
of 2D for SLG (solid circles) and 2D main components for 
BLG (open circles) as a function of n. 



further enhancing the phonon hardening. It is also pos- 
sible to demonstrate that, for T and 5 — > 0, the slope of 
APos(G) dyn just above Ep = 71 is double than that 
just below. Thus, the kink in Fig. 3 is a direct mea- 
surement of the interlayer coupling strength from Raman 
spectroscopy. 

In SLG the intensity ratio of 2D and G, I(2D)/I(G), 
has a strong dependence on doping[24 . Fig. 5a plots 
I(2D)/I(G) as a function of doping. For BLG we take the 
highest amongst 2Di,4 and 2T>2A- The SLG dependence 
reproduces our previous results [21]. However, we find an 
almost constant ratio in BLG. Fig. 5b plots the doping de- 
pendence of Pos(2D) in SLG, and Pos(2D 1A ), Pos(2D 2j4 ) 
in BLG. To a first approximation, this is governed by 
lattice relaxation, which explains the overall stiffening 
for hole doping and softening for electron doping 24J . A 
quantitative understanding is yet to emerge, and beyond 
DFT many body effects need be considered. 

To conclude, we have simultaneously measured the be- 
havior of optical phonons in single and bilayer graphene 
as a function of doping. In the latter, the G peak renor- 
malizes as the Fermi energy moves from the 1st to the 2nd 
subband, allowing a direct measurement of 71 ~0.4eV. 
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